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We present an infrared/optical study of the dynamics of the strongly correlated electron system 
YbIni_ a; Ag a ;Cu4 as a function of doping and temperature for x ranging from to 1, and T between 
20 and 300 K. This study reveals information about the unusual phase transition as well as the 
phases themselves. Scaling relations emerge from the data and are investigated in detail using a 
periodic Anderson model based calculation. We also provide a picture in which to view both the low 
and high-energy x -dependent features of the infrared data, including identification of high energy, 
temperature dependent features. 



When exploring the complex terrain of correlated elec- 
tron phases, one is often interested in the phase bound- 
aries separating systems which are controlled by dif- 
ferent physics. In the low doping regime (a;<0.2) of 
YbIn 1 _ :r Ag : ,.Cu4, a temperature-driven, first-order elec- 
tronic phase transition separates two phases which are 
characterized by widely disparate effective energy scales 
(Tr-) associated with the effects of hybridization between 
the localized Yb /-electrons and the itinerant conduc- 
tion electrons. This rapid change provides both an op- 
portunity to identify and isolate features associated with 
many-body hybridization physics and raises fundamental 
questions regarding the origin of the transition itself. 

The valence transition in YbInCu4 has been examined 
with a variety of experimental techniques 1-2-3 4*SiSii, and 
exhibits some similarities^*^ - to phase transitions ob- 
served in certain rare-earth intermetallic compoundsii*i£. 
Samples of high quality can be made, and provide the 
only known example of a valence transition occurring at 
ambient pressure in a stoichiometric sample. The synthe- 
sis of high quality single crystal samples has allowed the 
detailed study of the phase diagram of YbIni_ x Ag x Cu4, 
and continued effort has revealed a tendency of the sys- 
tem to order ferromagnetically as the phase transition 
temperature is driven to zero either by pressur e) 13 ! 14 or 
Y dopingi 3 * 1 ^. The presence of interesting phase com- 
petition in a system where the sample chemistry is well 
controlled enhances our interest in this correlated elec- 
tron system. 

Figure^shows a schematic representation of the phase 
diagram of Yblni-zAg^C^. At low doping, a line of 
first order phase transitions separates a low tempera- 
ture, mixed-valent phase and a high temperature lo- 
cal moment phase. In YbInCu4 (x=0) at high tem- 
perature (T > T v — 42 A"), the magnetic response ex- 
hibits a Curie- Weiss form with magnitude appropriate to 
j = 7/2 moment of Yb and a small Weiss temperature 16 
®w ~ —13 if, which constrains the effective Kondo 
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FIG. 1: Schematic phase diagram of Yblni-^Ag^Cu^ 
Lighter shading indicates a larger Tk- 



temperature appropriate in this high temperature range 
to be comparably small in magnitude. The resistiv- 
ity in this temperature ranged is very high for a metal 
(p^ c ~ 150/iO-cm) and together with Hall measurements 
indicate a very low (hole) carrier concentration^. Lower- 
ing the temperature through the first order phase tran- 
sition has dramatic effects on both the spin and charge 
response. For T < T v , the magnetic susceptibility drops 
considerably into a temperature independent (Pauli para- 
magnetic) form, with a magnitude indicating a dramatic 
increase in the Kondo temperature. In concert the car- 
rier concentration increases markedly and the resistance 
drops 2 . Extensive experimentation has revealed substan- 
tial changes in specific heat^i^*^, elastic constant a 6 ' 17 , 
single-particle (photoemission)i2i22i 2 i22i 2 £, two-particle 
(optical )2 4 *2Si2Si2 7 . and neutron^ spectroscopies. 

Changes associated with the phase transition in 
YbInCu4 are generally interpreted as an effective screen- 
ing of the localized Yb 4/ moments in the low tempera- 
ture phase. This is associated with an order of magnitude 
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increase in Kondo coupling scale (Tk-~17 K for T > T v ; 
T K ~300K for T < T v ). Universal scaling observed 4 -^ 
in the B-T plane was shown2___i to be consistent with 
the scenario of a transition-induced screening of the Yb 
moment. 

Doping Ag on the In site serves to stabilize the low 
temperature phase, driving the phase transition temper- 
ature upward—. Experimentally, the transition becomes 
less sharp with doping and it is estimated 2 that the tran- 
sition disappears completely at the critical concentration 
x c ~ 0.2. Concentration values x > x c display behav- 
ior more typical of heavy fermion systems, where co- 
herent heavy metallic quasiparticles responsible for the 
low energy physics have an enhanced mass as high as 
m* ~ 50m e . 

Much of the phenomenology of YbI1ii_j.Agj.Cu4 is re- 
flective of periodic Anderson model (PAM) physics, how- 
ever, the phase transitions known to exis1j2_in the context 
of the PAM are associated with a competition between 
magnetic order and Fermi liquid phases and occur in a re- 
gion of weak coupling. This region is far from the mixed- 
valent fixed point, where the coupling is strong or inter- 
mediate, and these transitions therefore seem unlikely to 
be directly relevant to YbIni_ I Ag x Cu4. In the region of 
intermediate to strong coupling the PAM exhibits grad- 
ual crossover to a coherent state as temperature is re- 
duced, and does not, to our knowledge, include phase 
transitions. One may therefore ask whether the line of 
phase transitions occurring at low x in YbIni_ :r Ag x Cu4, 
which are electronic and not structural in origin, im- 
ply the need for an additional interaction term in the 
minimal model Hamiltonian beyond those normally in- 
cluded in the PAM. An intriguing possibility is that 
electron-electron interactions, which can be very effective 
in bringing about transitions in Mott-Hubbard systems, 
may play a role in YbIni_ 2 .Ag x Cu4. Investigations along 
these lines have been pursued by Giamarchi et and 
Zlatic and Freericka 33 ! 34 ' 35 . 

In this paper, we present infrared-frequency dynam- 
ics derived from measurements of the reflectivity of sin- 
gle crystal YbIni_ x Ag :E Cu4, with Ag compositions x=0, 
0.3, 0.5, 0.75, and 1.0. This systematic study allows 
us to identify and explore scaling behavior of the in- 
frared signature of the Kondo resonance. In the x—0 sys- 
tem (YbInCu4), these measurements^ revealed a mid- 
infrared excitation prominent primarily in the low tem- 
perature phase212_ (T < T v = 42 K). This may be un- 
derstood as a consequence of the presence of hybridized 
quasiparticle bands2 4 *2____£i. Here we explore in detail 
the phenomenology of this excitation. Within the low 
temperature phase we use modeling to investigate the 
scaling properties of this infrared excitation. In addi- 
tion, we present interpretations of the higher frequency 
features based on systematic trends of spectral features 
with Ag substitution (hole doping) and propose an en- 
ergy level scheme in which to view the high-frequency 
spectroscopic results. 

Our measurements cover the frequency range from 



40 cm -1 to 50, 000 cm -1 with detailed temperature de- 
pendent data taken between 40 cm -1 and 23, 000 cm -1 . 
Each sample was mounted into a recess custom-machined 
into a brass disk. Each disk was polished (minimally) 
such that a smooth, flat sample surface sits parallel to 
a similarly Hat, polished portion of the brass mounting 
disk. Ag film was then evaporated on the bare brass 
surface, creating a reference mirror for use in the tem- 
perature dependent reflectivity measurements. Soon af- 
ter evaporation, the sample disk was mounted inside a 
continuous-How He 4 cryostat with a custom switching 
mechanism designed to expose either Ag or sample sur- 
faces to the beam. To insure that the reference mirror 
and sample surface are coplanar and Hat, a laser align- 
ment procedure was used on the mounted sample disk. 
Frequency dependent spectra were then taken using a 
combination of Fourier transform and grating spectrome- 
ters. At each temperature and composition, spectra were 
taken for Ag and sample exposed through the cryostat 
window and the the ratio at each frequency used to de- 
termine the reflectivity R(lu). 

In addition to the temperature dependence, room tem- 
perature reflectivity spectra were also taken in the range 
12,500 cm -1 < ui < 50,000 cm -1 with an optical setup 
which measures absolute refelctivity without the use of 
a reference mirror. The spectra taken by this method 
and those spectra taken using the cryostat (and with Ag 
reference) generally agree well in the region of overlap 
(12,500cm -1 < uj < 23,000cm -1 ). Through this com- 
parison, we estimate the uncertainty of the absolute value 
of the reflectivity in the temperature dependent measure- 
ments to be less than 1%. 

A Kramers-Kronig transform is applied to the mea- 
sured reflectivity in order to determine the frequency de- 
pendent reflection phase shiffe2____i. The magnitude and 
phase of the reflectivity are then used in order to deter- 
mine the dynamical conductivity u\{uj) (and the dielec- 
tric function ei(u)). 

For the purposes of the transform, Hagen-Rubens ter- 
minations (1 — R(u>) oc y/uj) are used below 40 cm -1 . 
At high frequency (above 50, 000 cm -1 ) each reflectiv- 
ity spectrum is extrapolated to a common value of 
0.08 at 120,000 cm -1 and then continued as a con- 
stant to 200,000 cm -1 . Between 200,000 cm -1 and 
400,000cm -1 , an w -2 form is used for R(w) to repre- 
sent the non-constant reflectivity expected as a result of 
deep core level excitations^. At still higher frequnecies, 
the free-electron form ui~ is assumed. 

We have experimented with a number of termina- 
tion protocols including other common values and co- 
alescence frequencies as well as constant extrapolations 
above 50,000 cm -1 . These show convincingly that our 
results regarding trends in the x and Tk dependence of 
(Ji(u>) below 10, 000 cm -1 are not significantly influenced 
by any of the extrapolations above 50,000 cm -1 . The 
conductivity above 10,000 cm -1 but below 20,000 cm -1 
is influenced by the detailed extrapolation by about 5% 
for reasonable extrapolation protocols. 
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In YbInCu4, where a first-order phase transition oc- 
curs at finite temperature, we have been careful to con- 
trol the potential hysteresis effects 2 ^ usually associated 
with first order phase transitions by keeping the ther- 
mal cycling through the phase transition to a minimum 
while collecting the optical data. We found that hys- 
teresis effects are observable in the infrared spectrum, 
and are especially pronounced in the frequency range 
1000 < lj < 7000 cm -1 (see Figure |2J- For this reason, 
our data for several ranges were retaken on fresh samples. 



I. RESULTS 
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FIG. 2: (Color online) The real part of the optical conductiv- 
ity of YbInCu4 at various temperatures. 

Figure shows the frequency-dependent reflectivity 
and infrared conductivity <t\(uj) of YbInCu4 at temper- 
atures below and above the T v ~ 42if phase transi- 
tion temperature. At high temperature (250 if), ct±(cj) 
consists of a narrow free-carrier (Drude-like) contribu- 
tion clearly seen in the far-infrared spectrum, due to 
the presence of mobile carriers. This interpretation is 
consistent with HalsS and resistivity measurements 2 . 
uj ~ 6000cm -1 marks a clear onset of a set of strong 
interband transitions extending upward into the visible 
spectral range. 

At the lowest temperature (20 K), a significant de- 
crease in conductivity in the interband region {ui > 



4000 cm -1 ) accompanies a substantial increase in the 
conductivity below 4000 cm -1 in the form of a well- 
defined peak centered around 2000 cm -1 . The develop- 
ment of this peak is highly correlated with the phase 
transition (7^=42 if) and is a prominent and essential 
feature of the low-T phase of YbInCu4 . The connection 
between this peak and the physics of hybridization is an 
important theme of the present work. 
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FIG. 3: (Color online) The real part of the optical conductiv- 
ity of YbIn.7Ag.3Cu4 at various temperatures. 

Figure [3] shows optical data for the a;=0.3 system. 
Trends similar to those in low temperature YbInCu4 are 
apparent, with a depletion of weight around 8000 cm -1 , 
accommodated by a replenishing at lower frequencies. 
Vestiges of the 2000cm" 1 peak persist at temperatures 
as high as 300 if, an energy scale much higher than the 
extrapolated value of the phase transition temperature 
for this composition T„(~100if ), but still lower than the 
Kondo scale appropriate to the low temperature phase 
T K ~ 360if. 

Figure 0] shows the similar plots for YbAgCu4 . The 
behavior displayed by the conductivity shows behavior 
more typical of a heavy fermion system, with moder- 
ate temperature dependence and spectral weight which 
is approximately conserved below 0.5 eV. There is a con- 
tinuous evolution of the 2000 cm -1 feature, activated by 
crossing the valence transition at low x, and the more 
familiar phenomenology found in the low frequency re- 
sponse of the heavy fermi system YbAgCu4. 
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FIG. 4: (Color online) The real part of the optical conductiv- 
ity of YbAgCru at various temperatures. 



FigureElprofiles the x dependence of oi(ui) at low tem- 
peratures. The 2,000 cm -1 peak in YbInCu4 undergoes a 
complex shifting behavior as x is increased, blueshifting 
slightly when £=0.3, then redshifting upon further dop- 
ing, reaching a minimum peak frequency when x=0.75, 
before blueshifting again as x continues to 1. The 
strength of this feature is also influenced by a; in a non- 
trivial way, discussed further below. 

At higher frequency, the large hump feature centered 
on 11,000 cm -1 in YbInCu4 monotonically blueshifts and 
decreases in overall strength as x is increased. Further, 
inflection points around 6000 cm -1 and 9000 cm -1 red- 
shift slightly upon doping to the x=0.3 and cc=0.5 sys- 
tems. These inflection points are not discernable for 
a;=0.75, but reappear at low frequency (3000 cm -1 and 
5500 cm -1 ) in the x=\ system. 

The cc-dependent profiling of the low temperature con- 
ductivity is an important part of our experimental re- 
sults, allowing clear identification of systematic changes 
of the low temperature electrodynamics as a function of 
an external control parameter. The discussion begins 
with identification of systematic trends in the 2000 cm -1 
feature followed by a discussion of interband features in 
a later section. 
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FIG. 5: (Color online) Infrared conductivity ai(uj) versus ui 
at 20 K for the five x values studied. Also shown are the 
components of a Lorentzian-Drude fit as described in the text. 
The dashed vertical lines are a guide to the eye and correspond 
to the center frequencies of the fit components when x=0. 



II. RESULTS, x DEPENDENCE 

In this section, we focus on the x dependent changes 
of the 2000 cm -1 feature in order to explore the signifi- 
cance of this behavior to the correlated electron physics 
of hybridization and the periodic Anderson model. First, 
we quantify the x dependent trends of Figure Eland draw 
important conclusions through the comparison to previ- 
ously published thermodynamic data. The inferred rela- 
tionships are then explored below in sections llVI and iLLTl 
where we consider an interpretation based on the periodic 
Anderson model. 
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FIG. 6: (Color online) (a) The characteristic frequencies, u) v k 
(triangles) and Uth (boxes), and (b) the Kondo temperature 
as a function of x. (c) The spectral weights, n(4000 cm -1 ) 
(boxes) and n(6000 cm -1 ) (circles) and the results of a com- 
bined Drude/Lorentz fitting (triangles) of the low frequency 
conductivity, (d) The low temperature susceptibility (open 
circles) and Sommerfeld coefficient (solid circles) as a func- 
tion of X. Tk in (b) is related to x(0) in (d). Spectral weights 
are computed assuming a band mass of 4m e . 



FigureEK shows the frequency of the 2000 cm -1 feature 
versus x as determined in two ways. The black triangles 
mark the frequency of the peak in a\(ui) (also marked 
in Figure [SJ. Alternatively, a threshold frequency can 
be extracted from a fit of the conductivity to a calcula- 
tion based on the low energy dispersion of the periodic 
Anderson model (PAM), discussed below in Section llVl 

In addition to examining the frequency of the peak 
as a function of x, we can also look at the strength of 
the 2000 cm - 1 feature. We quantify this characteristic 
through the spectral weight, defined as the integrated 
intensity of <7i(w) over a low frequency interval: 



n(uj) = ( cri(o/)du/ 
ne* J + 



(1) 



where m represents a bare band mass. The lower limit 
is chosen to be nonzero (0 + = 50 cm -1 ) in order to 
exclude from the strength estimate the comparatively 
minute contribution of the free carrier (Drude) response. 
The upper limit of integration is chosen to encompass 
the 2000 cm -1 peak without including the cc-dependence 
of the high frequency interband contributions. Neither 
integration limit is critical; in fact a lower limit of and 
upper limits between anywhere between 3000cm -1 and 
8000 cm -1 produce similar x dependence. n(4000 cm -1 ) 
and n(6000 cm -1 ) are shown in Figure^];. 

As an alternative to this simple integral calculation of 
the strength, we can fit the complex conductivity (a = 



The constituents of the fits include: one narrow Drude 
(Dl, r ~ 10 — 40 cm -1 ) contribution to represent the free 
^carrier peak, a wide Drude (D2, T > 800cm -1 ), and a 
,'gLorentz oscillator (LI) in the vicinity of the 2000 cm -1 

^feature which, as we discuss below, relate to the Kondo 
3 ... 
j ^resonance. There are in addition two Lorentz oscilla- 
tors (L2 and L3, around 7,300 cm -1 and 11,000 cm -1 for 
x=0) to represent the infrared interband conductivity, 
and two wide Lorentz oscillators at ultraviolet frequen- 
cies (u> > 30,000 cm -1 ) representing the conductivity in 
that range. These fit components are labeled in Figure 

Previous work^i has explicitly demonstrated that the 
calculated lineshape of the optical signature of the Kondo 
resonance is intrinsically non-Lorentzian, and further- 
more demonstrated the viability of fits which combine 
Drude and Lorentz terms to represent the Kondo res- 
onance. The combination of the contributions D2 and 
LI reasonably fit the conductivity in the range of the 
2000 cm -1 peak. The strength from that combination 
is shown by the triangles in Figure Efc. This determi- 
nation of the strength exhibits an x dependence similar 
to the simpler integral representations of the strength. 
This makes one confident that the n versus x dependence 
shown here is an essential characteristic of the data, and 
independent of any of the detailed choices we have made 
in the analysis, n determined by the methods discussed 
here is presented in Figure EJ:. 

It is interesting to compare the ^-dependent trends in- 
ferred from previously published thermodynamic mea- 
surements with those from our optical data. Fig- 
ure EJd shows the Kondo temperature of low-T phase 
YbIni_ 2; Ag a ;Cu4 as deduced by Cornelius 3 - et al from fit- 
ting the measured magnetic susceptibility to the numer- 
ically calculated result of the j — | Coqblin-Schrciffcr 
model 3 ^. Figure Eli shows that the trend in the x de- 
pendence of x(T = 0)S generally agrees with the corre- 
sponding trend in the Sommerfeld coefficient 7, implying 
a Wilson ratio within 10% of the value (1Z = |) expected 
for a j = I Anderson impurity 3 ^ 3 -. Thus the Tk values 
inferred from the susceptibility analysis reasonably rep- 
resent the effective energy scale relevant to the onset of 
strong-coupling Kondo physics within the low tempera- 
ture phase. 

The complicated x dependence of Tk is not under- 
stood, and may be the result of an interplay of band 
structure, chemical pressure, screening, disorder, and 
other many-body effects. While resolving the detailed 
cause of this complex x dependence presents a subject for 
future work, out focus here is the relationship between 
the x dependence of Tk and that of the optical data; the 
similar form of the x dependent electrodynamic (Figures 



6 







— — -% 






2V 

/ F~ 
/ 

f 









lr<l»ici n~ 
Kp 



FIG. 7: (Color online) The PAM dispersion relations. Vertical 
arrows indicate possible optical transitions. The horizontal 
dashed lines represent Ef and e/. The light diagonal line 
represents the unrenormalized dispersion of the conduction 
carriers. 

E^, and Et) and thermodynamic quantities (Figures |Hp 
andHJl) indicates a common source. We show below that 
this is rooted in the strongly correlated electron physics 
of hybridization. 



III. MODELING, PAM DISPERSION 

We can make progress toward eliciting the relation- 
ships suggested in Figures by examining an interpreta- 
tion of the 2000 cm -1 feature in the context of the pe- 
riodic Anderson model (PAM). Complementary to the 
work of other author which focus the rig- 
orous techniques of many-body theory directly toward 
the underlying Hamiltonian, we will use a simplified ap- 
proach based on effective low energy (near- Ep) PAM dis- 
persion relationship. These dispersion relations (e + and 
e~ in Figure provide the basis for a simple unifying 
picture in which much of the low energy phenomenology 
of heavy fermion materials can be viewed, including the 
mass enhancement, aspects of magnetism, and transport 
measurements^^. In Figure the vertical extent of 
the plot is of order 1 eV and the (singly occupied) /- 
electron level is below the bottom plot boundary. The 
light dashed lines indicate the bare (unhybridized) con- 
duction electron dispersion. 

In a system with no hybridization, the conduction 
electrons are the dominant influence on the transport 
properties such as thermopower and resistivity, while at 
the same time provide a temperature-independent Pauli- 
paramagnetic contribution to the magnetic susceptibility. 
The /-electrons, on the other hand, are localized and as 
a result contribute very little to the transport proper- 
ties, but play a major role in magnetism, contributing a 
Curie 1 /T term to the susceptibility. Inclusion of the hy- 
bridization and on-site Coulomb repulsion terms compli- 
cates this independent particle picture considerably and 
the new eigenstates become nontrivial admixtures of the 



states of pure / and conduction electron character. 

The PAM quasiparticle dispersion relations provide a 
venue through which to explore the effects of nonzero 
hybridization on the phase space of excitations. At ener- 
gies far from the chemical potential, the upper and lower 
bands, e + and e~, follow closely the unrenormalized free 
carrier dispersion. There the main effect of the inter- 
action and hybridization is to provide a channel for re- 
laxation of the conduction states, i.e. a broadening of 
the spectral function along the dispersion curves. At 
lower energies, these relaxation effects are reduced by 
the phase-space-constraining presence of a filled Fermi 
sea, however, renormalization also opens the Fermi sur- 
face (kp are ^ — > kp and the bands flatten) to accommo- 
date the /-electron weight projected up to the Fermi 
level. This reorganization of the bands in the vicinity 
of the Fermi level is due to many-body interactions, the 
strength of which is characterized by the parameter V. 
The resultant narrow peak in the density of states, called 
the Kondo, or Abrikosov-Suhl, resonance^ is central to 
the understanding of heavy fermion and mixed-valent 
phenomenologjs^ihi^. 

From the point of view of optical probes, a key effect of 
renormalization on the optical response is to create the 
possibility for vertical transitions from filled states below 
Ep, across a direct gap, and into unoccupied levels above 
Ep , as illustrated by the vertical arrows of Figure[7] The 
threshold for these transitions (short arrow) occurs at a 
frequency u) — 2V, where V is the hybridization strength 
renormalized by the on-site /-electron repulsion. This 
energy scales with the Kondo temperature a£ 48 - 51 - 52 ' 53 - 54 

V = ^JTkB (3) 

where B is a model parameter related to the conduction 
electron bandwidth. At threshold, the nesting condition 
for the upper and lower bands is met (Vj^e" 1 " = Vj^e - ), 
leading to a very high joint density of states for vertical 
quasiparticle transitions and hence a strong peak in the 
conductivity. 

At frequencies larger than the threshold frequency 
there are two distinct contributions to the conductiv- 
ity: one originating from levels inside the unrenormal- 
ized Fermi surface {k < below); the other from the states 
occupied as a result of renormalization, i.e., outside the 
unrenormalized Fermi surface (fc> below). An example 
of two such transitions with the same frequency uj are 
indicated by dashed arrows in Figure Transitions in- 
volving both of these sets of quasiparticles are important 
and must be counted independently in the determination 
of the total optical conductivity, as discussed further be- 
low. 

In addition to the threshold frequency, another, higher 
frequency scale appears which is relevant to the elec- 
trodynamic response. This higher frequency scale cor- 
responds to the vertical transition (long arrow, Figured) 
which occurs from states on the Fermi surface (i.e. the 
locus of points which divides the set of occupied and un- 
occupied k states). Vertical transitions involving higher 
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k states cannot occur because both initial and final states 
are unoccupied when k > hp, and hence one expects a 
drop in the conductivity at this frequency. For a lin- 
early dispersing conduction band, direct calculation re- 
veals that the frequency of the last allowed transition is 
equal to 



FS 



- 2 



(4) 



The identification of if with Tk (discussed further be- 
low), together with Equation |31 implies that this scale is 
of order the conduction electron bandwidth, B. For high 
energy transitions, band edge final states can be reached 
and the linear approximation to the conduction band is 
likely to become poor. Quasiparticle transitions in this 
frequency range may be influenced by the details of the 
underlying band structure. 



IV. MODELING, KUBO-GREENWOOD 
ANALYSIS 



With these considerations of the phase space for opti- 
cal transitions in mind, we can proceed with an analysis 
using the Kubo-Greenwood formula 38 : 



CTi (to) 



ire 



^JDOS u ,(uj)\p iti ,\ 



(5) 



where \p e ^> \ denotes the dipole matrix element connect- 
ing electronic bands £ and £', and JDOSe^'(u>) is the 
corresponding joint density of states. Applying this for- 
mula to hybridizing quasiparticles (as though they were 
electrons) allows an exploration of the phenomena of the 
mid-infrared conductivity in the context of the PAM. In 
that case, the two relevant bands are e + and e — , which 
in this model approach leads to: 



CTpam(^) 



dS 



4ir 2 m 2 u J Ae=w |V k (e+ - e 



(6) 



where |p_| | is the matrix element for transitions be- 
tween the two bands under consideration. 

In the case of a spherical Fermi surface, the integrand 
is constant and equation JSJ) simplifies to 



y pain 



(CO) 



47Tfc 2 | P+ ,_| 



47T 2 m 2 uj ^— ' 

k' — k^ ,fc> 



\d k (e^ 



(7) 



k=k' 



It is useful to consider the case of constant optical matrix 
elements |p + _|, meaning that every filled k state transits 
to the corresponding upper band state with equal proba- 
bility. One can now obtain a model lincshape from Equa- 
tion using the explicit PAM dispersion relations^ 3 ^ 8 *^ 



E F + e~ f + 6 k ± yj (Ep + if ~e k y+ W 2 
e 7. 1 ( 8 ) 



where Ep is the Fermi level, V is the renormalized hy- 
bridization strength and if is the /-level position renor- 
malized by on-site /-electron repulsion. This latter quan- 
tity defines the scale of the low-energy physics and is com- 
monly identified with the impurity Kondo temperature, 
Tk (discussed further below). 

Putting Equation |H1 into the conductivity takes the 
form: 



for 2V < to < 



2e 2 [p +| _l 2 (E F + i f ) 2 + to 2 - W 2 
m 2 \Ve kF \ 3 ir 



Vu 2 



W 2 



(9) 



and 



Cpam (k-0 



e 2 |p+,J 2 (Ep + if- Vlo 2 -W 2 ) 2 



TO 2 |Ve feF | 3 7T 



y/uo 2 - 4V 2 



(10) 



for to > 



Figure [SJi shows this lineshape for two 

sets of V and if values. 

The rather idealized lineshape generated by these con- 
siderations is extremely sharp and a meaningful compar- 
ison with the data requires addressing the effects of re- 
laxation, which were thus far neglected in our treatment. 
To this end, we convolute this idealized lineshape with 
a Lorentzian function, keeping the half width A as an 
adjustable parameter when performing fits to the mea- 
sured conductivity. Examples of fits produced using this 
procedure are shown in Figure [Bp. 

The broadening parameter A addresses the finite width 
of the spectral function along the dispersion curves of 
Figure and therefore provides a measure of the statis- 
tical time over which a typical quasiparticle decays. The 
numerical values for A obtained from our fits range from 
0.12 eV to 0.17ey, with an associated time scale for de- 
cay in the range r = h/ A=5Aps to 3.8ps, respectively. 
These lifetimes estimates are in good agreement with the 
quasiparticle lifetime of members of this class of mate- 
rials (YbAgCui), as measured directly by Demsar et al 
in pump-probe experiments^ of electron-hole relaxation 
lifetime. The extraction of this parameter from the con- 
ductivity data is a meaningful consistency check on the 
method developed here. 

We now take a moment to consider the possible in- 
fluence of k-dependent matrix elements associated with 
the composite nature of the hybridized quasiparticles on 
the electrodynamic response. The hybridizing quasipar- 
ticles are composite admixtures of excitations with both 
/ and conduction electron character. The regions of the 
dispersion which are flatter correspond to quasiparticles 
with a large amplitude of / admixture whereas regions 
which follow more closely the bare conduction disper- 
sion are dominated by conduction character. Generally 
speaking transitions among the bare states are not all 
equally probable (i.e. |p cc | ^ \P c f \ IP//I)- Thus one 
expects that the optical transition rate for quasiparticles 
may exhibit some dependence on k, which goes beyond 
our earlier assumptions. 
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FIG. 8: (Color online) (a) The idealized conductivity (Equa- 
tions [0] and Hot for V = 93 meV and e/ = 9 meV (solid) 
and = 40meV and e/ = 2.5 meV (dashed). Inset con- 
trasts the lineshapes derived with constant matrix elements 

p H | ( Equations l9l and II OH . and those with the coherence 

factors discussed in the text, (b) The same idealized con- 
ductivity curves as (a), Lorenztian broadened with widths 
Ao = 0.16 eV and A0.75 = 0.125 eV for comparison to 
the measured conductivity of YbIni- I Ag a; Cu4 with a;=0 and 
2=0.75, respectively. 
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FIG. 9: (Color online) Plot of uj p k and uith versus %/Tk". The 
dotted lines represent Eq. Q and the slopes are 1.6 eV? and 
l.leVi. 



conductivity at high frequency should fall to zero much 
more quickly than in the constant matrix element case^. 
The line shapes with and without coherence factors are 
contrasted in the inset of Figure 00 

This approach may be too nuanced because in a gen- 
eral mixed valent system, bare /-to-conduction state 
transitions can be appreciable, and in fact are expected 
to be important in Yblni-^Ag^Cru. This is because the 
conduction band states are derived primarily from Cu-In- 
Ag p and d orbitals^, whereas the / electrons sit on Yb 
sites. This physical displacement between the underly- 
ing orbital states is manifest in the banded states through 
nonvanishing dipole transition matrix elements, |pj c |££. 
Thus, in YbIni_ a; Ag a ;Cu4, transitions involving the flat- 
ter portions of the quasiparticle dispersion, which are 
dominated by /-like character, can be expected to pro- 
vide a considerable contribution to the optical strength, 
thus we feel that Ea uat ions 1^1 and HIH are more applicable 
to Yblni^Ag^Cuzt. 



One can obtain a relatively simple model with k de- 
pendence by assuming transitions among bare states only 
occur between conduction electron initial and final states 

(|Pccl 7^ 0)|P c /l = IP//I = 0)' an d using that to cal- 
culate transition rates between the hybridized quasipar- 
ticle bands. To model the non-constant admixture of 
states, we use the coherence factors of the resonant level 
model^SS (aka Fano- Anderson^), 

1 , . . 

u \c „ = — =, v\, „ = — 1 (11) 

The approach follows as before however the JDOS in- 
tegral (Eqn. H3) picks up a factor u \ La v ] iiy = V 2 /uj 2 as- 
sociated with these coherence factors. One thus obtains 
a model conductivity similar to Eg uat ions I§1 and 1 1 Ul but 
with |p + _| 2 replaced by \p cc \ 2 V 2 /ui 2 . The most signifi- 
cant effect of the inclusion of coherence factors is that the 



V. DISCUSSION, T K SCALING 

We now analyze the x dependent frequency of the 
2000 cm -1 feature (Figure EK), and its relationship to 
the Kondo temperature Tg (Figure Hip). Figure [5] shows 
tOpk, the maximum of the conductivity, and the threshold 
frequency uith, determined from the fit described above, 
plotted versus the square root of the Kondo tempera- 
ture. The complex x dependence of Lu t h and Lu p k (Fig- 
ure EJ) simplifies considerably when we plot these quanti- 
ties as a function of Tjc (Figure |5J). The emergence of a 
functional relationship between these quantities implies 
that the frequency of the 2000 cm -1 peak is controlled 
by the same physics that underlies the thermodynamic 
behavior. The square root dependence is evidence that 
hybridization physics plays a dominant role. 

With that in mind, the modeling developed in the 
previous section can be used to extract an estimate for 
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FIG. 10: Scaling relations for the strength of the 2000 cm - 
feature with Tk- Dark circles represent n(4000cm~ 1 ) and 
open circles represent ra(6000 cm - ) (a) shows the result of 
fitting the measured dependence with Equation 1141 and (b) 
shows the same data fit with Equation 1181 



the band parameter B (Equation |3J) by associating the 
threshold frequency, ui t h, with its PAM value, 2V. The 
slope of the line through the io th values in Figure El to- 
gether with Equation^ directly gives B = 0.30 eV. This 
value of B reflects the rate at which the frequency in- 
creases with Tk- This can be compared with expecta- 
tions based on density of states, as well as the rate at 
which the strength decreases with Tk, as discussed be- 
low. 



Within the approach developed in Section HVl the rela- 
tionship between Tk and the strength of the 2000 cm -1 
feature (n) can be addressed. We can obtain a closed- 
form result from our model calculation if we set the two 
sphere areas in equal to Ankp. This approximation 61 
avoids the effects of bare band structure details while 
including the influence of the strong cusp at 2V, which 
originates from the many-body physics of the PAM. Solv- 
ing for the wavevectors k > and fc< using the condition 
e + — e~ = to, and substituting the result into we ob- 
tain a model strength: 



2m 

vpam n 



2 V 



' pam 



(u))du> 



(12) 



-(£)• 



(13) 



This is essentially the area under the curves of Figure 
[S^,. Using Equation [3] and introducing the parameter c, 
defined by Tk = ce), we can express n pQm in terms of 
T K - 



Tlpam 




(14) 



Figure HOfa shows a least squares fit of this logarithmic 
scaling relationship to the data. 

The factor c relates the renormalized / level position at 
low energies and the Kondo temperature, Tk- The value 
of c is unambiguous in the Fermi liquid theory of the 
N(= 2j + l)-fold degenerate Anderson impurity model 
where, 



T K N 2 sin(7r/AT) cos(n/N) 
= k B T L — — 

C TT(N — 1) 



T K = w n T l . 
wn is the generalized Wilson number^ given by 



w N 



2nT(l + i; 



(15) 



(16) 



(17) 



For the j = 5 moment of Yb, N — 8 and c ~ 0.66. 

Using this value of c ~ 0.66 and fitting the measured 
result to Equation PHA gives B = 0.35 eV and 0.45 eV for 
n(4000 cm -1 ) and n(6000 cm -1 ) (Eq. 1), respectively. 
We consider these to be a reasonable agreement given 
the simplicity of our approach. These values of B, which 
are determined by the amount that n goes up when Tk 
goes down, are similar in size to our previous estimate of 
the parameter B — 0.30 eV determined from the amount 
that u>th goes down as Tk goes down. Therefore, in addi- 
tion to the agreement of experiment and theory regard- 
ing the direction of the Tk dependence of n and uJth, 
the sensitivity of the dependence of n and toth on Tk 
are in reasonable agreement with each other. In addi- 
tion, the numerical values for B are reasonable band- 
widths for the InAgCu ci-orbital derived band statesSSi^i 
of YbIni_ 2; Ag a ;Cu4. Thus we conclude that the observed 
dependences of both n and u>th on Tk is consistent with 
the predictions of the periodic Anderson model in mag- 
nitude as well as direction. 

An alternative strength estimate can be made using 
the coherence factor model of the conductivity, intro- 
duced in detail above. In that case, we replace |p + _| 2 
by \p cc \ 2 V 2 /u) 2 in EquationHJland repeat the steps which 
produced Equation 1141 giving a theoretical strength 



| Pcc | 2 fc F c 2 B - T K 
it 2 c 2 B + T K ' 



(18) 
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FIG. 11: (Color online) (a) The center frequencies of the tran- 
sitions L2 and L3 from (Figure 5), and (b) the strengths of 
the corresponding Lorentzian fit components are shown as a 
function of doping, x. A band-structure picture which we use 
to interpret the trends in (a) and (b) is shown for (c) YbInCu4 
and (d) YbAgCu4. The dashed lines illustrate k-conserving 
transitions which we associate with LI, L2 and L3, respec- 
tively. (LI is the Kondo resonance excitation.) 



Again using c = 0.66, fits to the data would yield B 
estimates 0.17 eV and 0.19eV, for n(4000cr7i _1 ) and 
7i(6000 cm -1 ), respectively. These values are in order-of- 
magnitude agreement with the band parameter estimates 
made above. This functional form is fit for comparison 
to the data in Figure ITUb. However, as discussed at the 
end of section IIV1 the /-to-conduction electron transi- 
tions are made allowed in this system and so the calcu- 
lation which includes these, represented by Equation 1131 
is probably more appropriate for Yblni^Ag^C^. 

VI. DISCUSSION, HIGH FREQUENCY 
INTERBAND TRANSITIONS 

We now consider the x-dependent trends in the high 
frequency (u > 6000 cm -1 ) conductivity. We will use 
the language introduced in Figure [S] (Section [H} asso- 
ciated with Lorentzian fitting of the conductivity data, 
focussing our attention on the features L2 and L3. In 
Figures ITTk and II lb we show the center frequency and 
strength of these high energy excitations as a function of 
x. This measured ^-dependence can inform our under- 



standing of the nature of the underlying states associated 
with these transitions and the density of states near Ep. 

In a textbook picture of metals and semiconductors, an 
important effect of doping is to add or remove electrons 
from a set of band states, thereby influencing the position 
of Ep. In YbIni_ ;E Ag 2 ;Cu4, increasing x from to 1 cor- 
responds to the net removal of 2 electrons (per formula 
unit) from the system, hence we expect the Fermi level 
to move downward in energy as x is increased. Relativis- 
tic band structure calculations addressing these changes 
have been carried out by Antonov et a&. In this sce- 
nario, optical features involving transitions from filled 
states just below Ep tend to weaken and move upward as 
x is increased and these states are emptied^. The dashed 
arrow of Figure ITTb illustrates one such transition, which 
we identify with L3. In this scenario, the width of the 
feature when x=0 suggests bandwidths of order ~leV, 
and the threshold for these transitions, approximately 
7000 cm -1 (^0.9 eV) when £=0, gives an indication of 
the overall energy position relative to Ep. Furthermore, 
the amount of shift with doping implies that the density 
of states in the neai-Ep region of the band structure is 
approximately 2/0.58 eV^S.be~~ eV' 1 /i.n. Both of these 
numbers are quite reasonable for conduction bands in 
rare earth and transition metal systems^. 

With this rigid band interpretation as a backdrop de- 
scribing the salient relationships between L3 and x, the 
correlated electron effects discussed in previous sections 
occur in addition. The renormalization discussed there 
dresses these bare states and the renormalized region 
of the band structure, including the Kondo resonance, 
tracks Ep, which moves downward with increasing x. 

We now turn to the systematics of the feature L2. This 
feature redshifts with x by an amount similar in magni- 
tude to the shift of L3 (Figure lllb ) . However, the nar- 
rowness and nearly x-independent strength of L2 does 
not lend as easily to a simple band interpretation. When 
considering an identification of the component L2, we 
point out the significant temperature dependence in the 
frequency region associated with L2, as shown before in 
Figures [3 El and El The temperature dependent interplay 
of spectral weight contained in the frequency intervals 
of LI and L2 naturally lead one to speculate that per- 
haps the same correlated electron physics controlling LI 
may also be relevant to L2. One is thus led toward the 
question of whether the presence of the feature L2 rep- 
resents a further phenomenon associated with hybridiza- 
tion physics. 

We noted in Section ILTT1 that in addition to the strong 
peak in the PAM conductivity associated with renor- 
malized band nesting, a second feature could appear at 
higher frequencies associated with the initial state en- 
ergy crossing Ep. In a linearly dispersing band model, 
this change occurs around the conduction electron band- 
width frequency (Equation In a more realistic band- 
structure, the filling fraction and band curvature details 
could influence the frequency and magnitude of this con- 
ductivity change. In particular, if the transition final 



11 



states are band edge states, then significant shifting with 
doping concentration could result. 

With these considerations, it is reasonable to suggest 
that L2 represents a conductivity feature arising from 
Fermi surface quasiparticles. The narrowness can then 
be attributed to the very long lifetime of the initial state 
(because it occurs on the Fermi surface), and weak x de- 
pendence of the strength arises because the energy struc- 
tures responsible for the associated transitions track the 
x dependent Fermi level, as opposed to becoming filled 
or depleted with x. The culmination of the identifica- 
tions suggested in this section are presented at a light 
and heavy doping concentration in Figures UTI and lllH . 
respectively. 

VII. DISCUSSION, THE PHASE TRANSITION 
OF YbInCu 4 

The physical mechanism of the phase transition in the 
lightly doped system remains elusive, however, progress 
has been made in understanding aspects of the phase 
transition. In previous sections, we have analyzed in de- 
tail the low temperature conductivity in terms of PAM 
renormalizations of few-band models. We now discuss 
the how this picture identifying optical transitions may 
be related to the interesting phenomenology displayed 
by YbIni_ a; Ag :E Cu4 including the phase transition at the 
low x. 

The Kondo volume collapse^ (KVC) model, which 
describes the complex interplay of the system volume, 
hybridization, Kondo temperature, and / level occupa- 
tion, seems adequate to describe the valence transition 
in elemental C o 12 ' 58 . This model, however, seems insuf- 
ficient to quantitatively describe the phase transition in 
YbInCu4, as evidenced mainly by the smallness of the 
volume change at the transition. Some author s^' ll-ty, 
have argued that a quasigap, or region of low density 
of states, exists in the bare band structure located just 
above the Fermi level in YbInCu4, and is important to 
the phase transition. This description is qualitative, but 
addresses the change in carrier density as well as the 
changes related to Kondo physics. In this approach, the 
presence of a quasigap makes the Kondo temperature 
very sensitive to the placement of the Fermi level and 
can help induce the phase transition in manner akin to 
the KVC model, but with the importance placed on the 
density of states dependence, rather than the volume. 

It is of interest to ask what consequences a quasigap 
scenario could have for the optical data in YbInCu4, 
where the conductivity in the range of 8000 cm -1 (L2 
in the fits) drastically displaces to lower frequency form- 



ing the Kondo resonance excitation (LI) discussed above. 
In the picture outlined in Figure ITTb . an upward shift in 
the Fermi level off of the bare band edge and into the 
quasigap could have drastic consequences for the optical 
response, possibly forcing the contribution L2 to shift to 
very high frequency, as the renormalized portion of the 
upper band is forced across the gap in response to a small 
shift in Ep . This possible interpretation of the temper- 
ature dependence indicates that further theoretical work 
directed toward investigating the generic physics of the 
PAM and the fate of the Kondo scenario in the context of 
rapidly varying band structure is needed to understand 
the complex temperature dependence of YbInCu4. 

Theoretical work on the Falicov-Kimball model2&2!i2£ 
attempts to make a quantitative connection with experi- 
ment using a particular many-body model which features 
a first-order phase transition. Aspects of this modeling 
seem promising, in particular the prediction of substan- 
tial temperature dependence of the high frequency op- 
tical conductivity, however this model does not yet ex- 
plicitly include the effects of hybridization. It is possible 
that a minimal model that describes the first order phase 
transition and also includes the Kondo physics may re- 
quire and extension of the periodic Anderson model to 
include Falicov-Kimball-type interaction. 



VIII. CONCLUSIONS 

Our results indicate that Tk scaling is present 
in the low temperature finite-frequency dynamics of 
YbIni_ a; Ag x Cu4 and can be addressed in the context 
of local-moment models. Furthermore, our data provide 
numerical estimates of key parameters necessary for the 
construction of a minimal theoretical model of the va- 
lence transition as well as pointing out salient features of 
the underlying bandstructure. Further work may be di- 
rected toward greater understanding of this low-T scaling 
behavior as well as unexplained temperature dependent 
behavior of lightly doped YbIni_a;Ag x Cu4. 
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